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ABSTRACT 

We present distribution functions and spatial correlations of the shapes of dark matter halos de- 
rived from Hubble Volume simulations of a ACDM universe. We measure both position and velocity 
shapes within spheres encompassing mean density 200 times the critical value, and calibrate small-N 
systematic errors using Poisson realizations of isothermal spheres and higher resolution simulations. 
For halos more massive than 3 x 10^** H^^Mq, the shape distribution function peaks at (minor/major, 
intermediate/major) axial ratios of (0.64, 0.76) in position, and is rounder in velocity, peaking at 
(0.72,0.82). Halo shapes are rounder at lower mass and/or redshift; the mean minor axis ratio in 
position follows (c/a)(M,z) = ci5,o[l - a\n{M/10^^ h-^MQ)]{l + with ci5,o = 0.631 ± 0.001, 

a = 0.023 ± 0.002 and e = 0.086 ± 0.004. Position and velocity principal axes are well aligned in 
direction, with median alignment angle 22°, and the axial ratios in these spaces are correlated in 
magnitude. We investigate mark correlations of halo pair orientations using two measures: a simple 
scalar product shows > 1% alignment extending to 30 Mpc while a filamentary statistic exhibits 
non-random alignment extending to scales ~200 Mpc, ten times the sample two-point correlation 
length and well into the regime of negative two-point correlation. Cluster shapes are unaffected by 
the large-scale environment; the shape distribution of supercluster members is consistent with that of 
the general population. 

Subject headings: cosmology: theory — large-scale structure of the universe — galaxies:clusters:general 



1. INTRODUCTION 

Clusters of galaxies signal the largest gravitationally 
bound dark matter halos in the universe. They are mildly 
aspherical systems that tend to be aligned by mergers 
directed by interconnecting filaments in the cosmic web. 

Since the early days of extragalactic astronomy, it has 
been apparent that clusters are generally elongated on 
the sky. Flatteni ng of these clusters due to rotation 
was ruled out (;Illingworth 1977^ and many assumed that 
the flattening was due to gravitational instabilities ex- 
pected in the top-down scenario of Zel'dovich (1978) and 
Doroshkevich et al. (1978). The groundbreaking work 
of Carter and Metcalfe (1980) showed that the aspher- 
ical shape of a cluster was connected with the velocity 
anisotropy of the orbits of cluster galaxies. Binney and 
Silk (1979) proposed that this anisotropy was due to tidal 
distortion from neighboring large-scale structure. 

Binggeli (1982) was the first to investigate alignments 

Electronic address: |skasun@umich.edu, evrard@umich.edu| 



of close cluster pairs. He studied 44 Abell clusters and 
found that galaxies separated by less than 30 Mpc show 
a strong alignment with each other and that the ori- 
entation of a cluster was dependent on the distribu- 
tion of surrounding clusters. Other observational claims 
of cluster alignments have been made, both with near- 
est neighbo rs and with other cluster s in the same su- 
nerrJuster Jwest 1989a[ Iwest 1989hl khee et al. IQoS 



Richstone et al. 1992HPlionis 1994^ . 



With some exceptions (for example, Struble and Pee- 
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confirms Binggeli's original results. West (1989a) uses 
48 Abell superclusters and finds a tendency for clus- 
ters within 60 h^^ Mpc to be aligned. Plionis (1994) 
measures alignments for 637 Lick clusters and finds 
strong alignments up to 15 Mpc with weaker align- 
ments out to 60 Mpc. West et al. (1995) finds a 
marked anisotropy for Einstein clusters extending out to 
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10 h-^ Mpc. 

Simulations have shown that spatial alignments, 
intrinsically pre dicted from the 'top-down' model 
i zei'dovich imh . are also seen in bottom-up 
(cold dark matter) scenarios, where halos form by 
mergers of smaller structures orga nized along fil- 



TABLE 1 
Model Parameters 



aments Istruble fc Ftaclas 1994 : 

West et al. 199,4 . Simulations also show that the 
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orientation of the major axis of a galaxy clus- 
ter is aligned with the direction of the last major 
event Jvan Haarlem fc van de Wevgaert 1993 



merger 



Sphnter et al. 1997|) . Since cluster alignments appear to 
be a generic outcome of gravitational instability, their 
use as a discrimin ant of cosmological mod els requires 
careful calibration l Onuora fc Thomas 200o() . 

In the halo model description of non-linear structure 
(Cooray and Sheth 2002), all matter is contained in 
bounded, spatially correlated regions (the halo popula- 
tion) that span a spectrum of sizes. Most instances of 
this model assume spherically symmetric halos, but more 
precise versions will need to take into account the spec- 
trum of halo shapes, including detailed internal structure 
(Jing & Suto 2002), as well the spatial correlation of the 
shapes of neighboring halos (Jing 2002; Faltenbacher et 
al. 2003). 

In this paper, we report measurements of shape statis- 
tics, including spatial (or mark) correlations of align- 
ments, derived from large samples of massive dark mat- 
ter halos extracted from Hubble Volume simulations. 
We investigate two flat-metric cosmologies (ACDM and 
rCDM) dominated by vacuum energy and dark matter, 
respectively. We focus on the former, more empirically 
satisfying model, but show results for the latter for com- 
parative purposes. Section 2 provides details of the sim- 
ulations and describes our method of finding principal 
axis orientations and magnitudes for mass-limited halo 
samples. In § 3, we use Poisson realizations of isother- 
mal spheres to estimate systematic error in mean shape 
measurement due to shot noise. We then present axial 
ratio distribution functions for mass-limited samples, and 
characterize the dependence of halo shape on mass and 
redshift. Section 4 presents mark correlations of ACDM 
cluster alignments and briefly examines the role of su- 
percluster membership on shape. A final section reviews 
our conclusions. 

Unless explicitly stated otherwise, the halo mass M 
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used throughout this paper is a critically-thresholded, 
spherical overdensity mass (Af2oo) expressed in units of 
1015 /i-iMq, with Ho/lOOkm s"i Mpc^^. 

2. CLUSTER SHAPES FROM SIMULATIONS 

We use two sets of simulations to investigate dark mat- 
ter halo shapes. The large Hubble Volume (HV) simula- 
tions provide statistical power but poor mass resolution, 
while higher resolution, but smaller volume, Virgo simu- 
lations are used for resolution tests. 

2.1. Simulations 

The HV simulations are a pair of gigaparticle N-body 
simulations cr eated using a parallel ve rsion of the Hydra 
N-body code i Macfarland et al. 199s[> . Random realiza- 
tions of two cosmologies with flat spatial metric are pro- 
duced: a ACDM model with = 0.3, fl\ = 0.7 and 
power spectrum normalization ag = 0.9 and a tCDM 
model with fl„i = 1, JIa — 0, and as — 0.6. The 
dark matter structure is resolved by particles of mass 
2.2 X 10^^ h~^MQ within periodic cubic volumes of length 
3000 and 2000 /i"! Mpc, respectively (see Table HI). We 
analyze z = and combined sky survey samples of galaxy 
cluster halos published in Evrard et al. (2002). The 
reader is referred to that paper for details of the simula- 
tions, including the process of sky survey creation. 

For resolution tests, we employ the 256'^-particle Virgo 
simulations of Jenkins et al. (1998). These simulations 
have order-of-magnitude improved mass resolution but 
much smaller samples than the HV simulations. 

2.2. Halo Finding Algorithm 

We define dark matter halos with a spherical overden- 
sity (SO) group finder that identifies as a halo the set of 
particles lying within a sphere of size r2oo , centered on a 
particle that represents a local density maximum filtered 
on a scale of 2 x 10^^ Ii^^Mq. The size measure r2oo is 
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the radius of the sphere within which the mean density 
is 200pc{z), with Pc{z) the critical density at redshift z. 
The total mass Af2oo lying within r2oo is the basic or- 
der parameter of the halo sample. As a result of ongoing 
merging activity, roughly 7.5% of halos are found to have 
overlapping r2oo spheres. We retain both members of an 
overlapping pair in the halo sample. 

2.3. Sky Surveys 

In addition to the traditional mode of fixed proper- 
time output, sky survey samples, consisting of data col- 
lected along the past light cone of hypothetical observers 
located within the computational volume, were also gen- 
erated by the HV simulations. We use two octant sur- 
veys (PO and NO) that cover tt/2 sterad and extend to 
z^a_^ = 1.46 for ACDM and 1.25 for rCDM, as well as two 
full-sky surveys (MS and VS) that reach Zmax = 0-57 for 
ACDM and 0.42 for tCDM. The octants sample struc- 
ture over the last 74% (ACDM) and 71% (rCDM) of the 
age of the universe, approximately a 10 Gyr look-back 
time. Halos in these surveys are defined using the SO 
method described above. 

2.4. Cluster Shapes 

Given the modest resolution of the HV data, we take a 
simple approach and estimate the shape of a halo using 
the moments of the material within r2oo ■ 

With respect to the center of the halo (defined by the 
local gravitational potential minimum), we compute the 
3x3 symmetric tensor Mjk 

Mjk = ^ ^XajXak (1) 

a 

where the j*'' component of the displacement vec- 

tor of particle a relative to the halo center and Nh is 
the number of (equal mass) particles in the halo. We 
diagonalize to find eigenvalues \j and unit eigenvectors. 
The sorted eigenvalues (Ai > A2 > A3) and vectors de- 
fine principle axes a, 6, c, with semi-major axis a = \/A]", 
intermediate h — •\/A2 , and c = VA3 of a triaxial ellipsoid 
that approximates the halo. We refer to the correspond- 
ing unit vector directions as a, b and c. 

Velocity moments are solved for in a similar way, us- 
ing a mean defined as the center of mass velocity of the 
material within r2oo- For compactness of notation, we 
define intermediate and minor axial ratios as follows 

b = b/a, (2) 
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Fig. 1. — Estimates of the effect of Poisson noise on measurement 
of minor (c) and intermediate (6) axis ratios. Halos must lie in the 
upper left portion of the plot, with spherical objects at (1, 1), oblate 
halos tending to b=l and prolate halos lying near the diagonal b = 
c. Contours show 68% confidence regions of axial ratios measured 
from 10,000 random realizations of isothermal spheres distorted 
to mean shape (0.65,0.80), shown by the cross, and resolved by 
(from outer to inner contours) 40, 160, 640, and 2560 particles. 
Histograms in the lower right give frequency distributions of the 
minor axis ratio c. As resolution degrades, the mean is biased to 
lower values and the dispersion grows. 

c = c/a. (3) 

We have confirmed that the halo population is oriented 
randomly with respect to the Cartesian coordinates of 
the simulation, as required by the cosmological principal. 

3. HALO SHAPE STATISTICS 

We begin this section by estimating the error on shape 
determinations due to shot noise. From Monte Carlo re- 
alizations of distorted isothermal spheres, we find that 
population mean values of c and b are accurate to a few 
percent if Nh ^ 100. We then present the joint dis- 
tribution function p(c, 6) at 2: = 0, and investigate joint 
position and velocity shape statistics. We follow with the 
mass and redshift dependence of the minor axis ratio. 

3.1. Resolution Tests 

To calibrate the error in mean shape due to the small 
numbers of particles used to resolve clusters in the HV, 
we carry out the following resolution test. Halo models 
are created using a set of Nh particles with an initially 
spherical, isothermal profile. To simulate the mean el- 
lipticity of the HV cluster samples discussed below, each 
cluster is then compressed along the (x, y)-axes by fixed 
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TABLE 2 

Resolution Test Values 



Nh (C) <Tg (b) 



amounts (0.65, 0.8). The particle moments are calcu- 
lated and used to estimate the shape of each cluster in 
the same manner as for the HV clusters (see section 2.4). 
Generating ensembles at several values of Nh, we cal- 
ibrate the bias in mean cluster shape as a function of 
mass resolution. 

The results for an ensemble of 10,000 clusters at each 
Nh are presented in Figure^ Mean axial ratio values (c) 
and (&) are presented in Table [3 for Nh ranging from 40 
to 2560 particles. In Figure ^ the measured shape dis- 
tribution function has a dispersion ranging from 0.09 at 
Nh=4:0 to 0.01 at Nh=2560, and the mean is offset from 
the input values by a bias that scales inversely with par- 
ticle number 0.024(iV;i/100)"^ We wiU use this calibra- 
tion to bias-correct the mass-dependent shapes presented 
in § 3.4. 

3.2. Shape Distributions 

To ensure a mean shape measurement biased by less 
than two percent, we impose a mass limit M200 > 
3 X 10^"* /i~^Mq, corresponding to 133 particles in the 
HV simulations. We note that Jing (2002) also finds that 
160 particles are sufficient for percent-level shape mea- 
surement. To keep the corrections small while maintain- 
ing a moderate sample size, the Virgo data are cut at a 
mass of 10^^ h~^MQ, equivalent to 1467 particles for the 
ACDM case and 440 particles for the tCDM case. The 
numbers of halos in each sample are listed in Table El 

The distribution of axial ratios ai z = for 82,967 
(ACDM) and 87,121 (rCDM) halos are given in Figure 
[3 The top row shows the distributions in position space 
while the bottom row shows velocity space. The location 
of the mean axial ratios ((c), (5)) is depicted by a cross 
for HV and an asterisk for Virgo data. 

Consistent with previous studies, we find galaxy clus- 
ter halos to be mostly prolate in shape, and somewhat 
rounder in velocity compared to position space. Table 



TABLE 3 
Halo Samples at z = 



Model Mmi„[h-^MQ] 

3 X IQi* 82967 

1014 353 

3 X 10^'' 87121 

10" 703 



Mass-limited samples, M200 > Af? 



ACDM tCDM 



J 1 , 1 , 1 1 1 1 1 1 l>+-Ht 
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0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 
c 

Fig. 2. — Contours of the joint probability distribution of axial 
ratios p{c, b) are shown with the minor-axis frequency distribution 
p(c) in the same format as Figure Left and right columns show 
the ACDM and rCDM models while upper and lower rows give 
position and velocity space distributions. Contour levels in the 
joint pdf are drawn at the 10th, 50th, and 95th percentiles of the 
enclosed distribution. Gaussian fits to the frequency distribution 
of c are shown for the HV simulations while histograms show c 
values derived from the Virgo simulations. Crosses and asterisks 
mark the locations ((c), (b)) of the mean shapes for the HV and 
Virgo models, respectively. 



El summarizes the z = shape data. The ACDM HV 
simulation has modal values {c,b)p^^f. — (0.64, 0.76) and 
(^Mpeak — 0.82) in position and velocity, while 

the rCDM model halos are more strongly ellipsoidal, 
with values (c,6);,,, = (0.59,0.72) and (5,6)^,,, = 
(0.70,0.80). 

The frequency distribution of minor axis ratio c, com- 
puted by integrating the joint pdf along the b axis, is well 
fit by a Gaussian for the HV samples. Distributions of c 
derived from the Virgo simulations, shown as histograms 



40 0.584 0.089 0.777 0.095 ACDM-HV 

160 0.635 0.051 0.794 0.061 ACDM-Virgo 

640 0.647 0.026 0.799 0.032 tCDM-HV 

2560 0.649 0.013 0.799 0.016 rCDM-Virgo 
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TABLE 4 
Mean Halo Shapes at z = ^ 















Model-Space 
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ACDM-P 


(0.635, 0.760) 


0.635 


0.086 


0.776 


0.096 


ACDM-V 


(0.719, 0.823) 


0.704 


0.080 


0.818 


0.088 


rCDM-P 


(0.593, 0.719) 


0.600 


0.087 


0.754 


0.102 


tCDM-V 


(0.698, 0.802) 


0.686 


0.077 


0.814 


0.086 



"Mass-limited samples of Table ISl 



in Figure|21 are generally in good agreement with the HV 
data. The distribution of ACDM minor axis ratios for po- 
sition is centered on 0.635 and has measured dispersion 
0.086. When corrected for Poisson error, the estimated 
mean is slightly larger (c) = 0.647 and the intrinsic dis- 
persion is estimated to be 0.076. While the statistical 
uncertainty in the mean is in the third significant digit, 
the level of systematic uncertainty is certainly larger. 

A crude estimate of systematic error is given by the 
0.012 Poisson bias correction for the mass-limited sam- 
ple. However, at the level of 0.01, a number of other 
effects on axial ratio come into play. Foremost among 
them is the detailed definition of a halo: its location, 
scale, and geometry. We use the common spherical over- 
density (SO) definition of halos as spherical regions cen- 
tered on local density peaks. Other viable approaches, 
such as friends-of-friends grouping or use of an ellipsoidal 
boundary can systematically shift shape measurements 
by many percent (Warren et al. 1992, Jing & Suto 2002). 
We leave it to future work to address such systematic ef- 
fects in detail. 

The finding that halos are rounder in velocity space 
may be at least partly due the effects of ongoing merg- 
ers. Mergers will scatter in velocity space first, followed 
by mixing and relaxation of particle positions. Another 
factor pushing in the same direction is that the gravita- 
tional potential that drives the velocity field is rounder 
than the density distribution. 

We do not attempt to formally fit the joint probability 
distribution p(c, 5), as doing so would requite a full de- 
convolution of the effects of shot noise. However, to give 
an indication of the joint pdf shape, we locate the peak 
in the intermediate axis conditional probability 



and record the modal intermediate axis ratio bmod as a 



TABLE 5 

Intermediate Axis Modal Ridge Line b^odic) 



Model Component r s 

ACDM Position 0.81 ± 0.02 0.24 ± 0.01 

Velocity 0.84 ± 0.03 0.20 ± 0.02 

rCDM Position 0.79 ± 0.02 0.25 ± 0.01 

Velocity 0.79 ± 0.03 0.25 ± 0.02 

Tit to 6„od(c) = rc + s. 




0.2 0.4 0.6 0.8 1 
|cos(0py) 



Fig. 3. — Cumulative distributions of the alignment angle be- 
tween the position and velocity major axes are shown for the 
ACDM HV clusters with M200 > 3 X lO^-* H-'^Mq (solid line) and 
ACDM Virgo clusters with M200 > 10^" /i-^M© (dashed hne). 
Half of the clusters in the ACDM HV simulations have an align- 
ment angle smaller than 22 degrees. 

function of minor axis ratio c. The lines in each panel of 
Figure 121 show fits bmod{c) — rc+s, with fit parameters 
given in Table [S] 

3.3. Position-Velocity Major Axis Alignment 

Since elongated orbits drive both position and veloc- 
ity anisotropy, a correlation between both position and 
velocity major axes is both expected and measured in 
simulations (Tormcn 1997). We quantify the alignment 
between position and velocity in a halo through the scalar 
product of its major axis eigenvectors 

cos(6lpy) = \a^-a^\. (5) 

Figure |21 shows cumulative probability functions of this 
statistic for the ACDM HV and Virgo models. 
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Fig. 4. — Joint probability density of the position and velocity 
minor axis ratios. Light contours show the 16th, 50th, 84th and 
99th percentiles of the enclosed distribution while dotted contours 
show the same for the ellipsoidal Gaussian fit described in the text. 
The bold contour highlights the median of the joint distribution. 



We find a strong alignment signal: half of all halos 
have position and velocity major axes aligned to better 
than 22 degrees for the HV ACDM model (21 degrees 
for rCDM). The higher resolution Virgo runs show even 
stronger position- velocity alignment, with median values 
of 15 degrees for ACDM and 14 degrees for the rCDM 
simulation. When the same analysis is carried out using 
only strongly ellipsoidal clusters (c < 0.635), the results 
do not significantly change. 

Our results are consistent with a higher resolution 
study of Tormen (1997), who finds from simulations of 9 
dark matter halos resolved by 20,000 particles that the 
position- velocity alignment angle is approximately 30 de- 
grees. 

To further define the relationship between position and 
velocity space, we show the joint distribution of minor 
axis ratios (c^,c^) in Figure 2| The likelihood is well 
fit by a Gaussian in an ellipsoidal coordinate r, where 
= (^)^ + (^)^andc± = ^^[{5^ ^{5P))±{5^ -{5^))] 
are principal component directions centered on the one- 
component means of the distribution. This fit is shown 
by the dotted lines in Figure^ The measured dispersions 
for ACDM are <t+ = 0.096 and cr_ = 0.069. The bold 
contour shows the median of the enclosed distribution. 

Since essentially all cluster detection methods use some 



power of projected mass density as a defining signal, the 
spatial- velocity alignments examined here will introduce 
an orientation bias in estimates of velocity dispersions 
for clusters lying close to the sample detection threshold. 
The magnitude of this effect will depend on the specific 
sample in question, but the general trend will be for the 
line-of-sight dispersion to overestimate its isotropic coun- 
terpart, by a fractional amount that could be ^ 10% for 
small samples of moderate signal-to-noise detections. Es- 
timating the effect for particular observational surveys is 
best done using Monte Carlo sample realizations, such as 
those performed for the 2MASS (Kochanek et al. 2003) 
and SDSS (Miller et al. 2004; Wechsler et al. 2004) sam- 
ples. 

3.4. Mass and Redshift Dependence 

Because mergers are directed along large-scale fila- 
ments in the cosmic web, the birth of dark matter ha- 
los is an inherently asymmetric process. As a merger 
evolves, dynamical relaxation will tend to drive a halo 
closer to isotropy. One then expects that dynamically 
younger clusters will be more strongly ellipsoidal than 
older ones. Equating dynamical age with elapsed time 
from a halo's formation epoch (defined, for example, us- 
ing the mass accretion history by Wechsler et al. 2002), 
we expect high mass halos at a given epoch will be more 
elongated than those of lower mass. Similarly, at fixed 
mass, high-z halos should be more ellipsoidal than their 
low-z counterparts of the same mass. 

With the large number of clusters in the z = sample, 
we first investigate the mass dependence at the present 
epoch. Figure [S] shows the dependence of mean minor 
axis ratio shapes in position (filled symbols) and veloc- 
ity (open symbols) on halo mass, using bins of width 
O.ldex. Error bars on the points give the uncertainty in 
the mean. Circles show the HV data while squares show 
the Virgo data. Small symbols show the raw mean shape 
measurements while large symbols correct for the bias in 
mean shape calibrated in section 3.1. There is generally 
good agreement between the bias-corrected mean shapes 
of the HV models and the values measured directly from 
the higher resolution Virgo runs. 

The mass dependence favors more ellipsoidal halos at 
higher mass. The mean minor axis ratios of both the 
position and velocity shapes vary by only ~ 5 percent 
as the mass changes by a factor of ten. For the ACDM 




Fig. 5. — The dependence of mean minor axis ratio on mass is shown for halo samples at z=0 of the a) ACDM and b) tCDM models. The 
upper curve with open symbols give velocity shapes while the lower curve with filled symbols show position. Small symbols show measures 
uncorrected for small-N effects while bigger symbols have Poisson corrections applied, as discussed in the text. Circles and squares are HV 
and Virgo model results, respectively, and the vertical dotted line marks the 3 X 10^'*/i~^Mq resolution limit of the HV models. Solid 
lines show the weak logarithmic dependence on mass, equation JSJ, with parameter values given in Table |S| Higher mass clusters are more 
strongly ellipsoidal than low mass clusters in both position and velocity. 



model, a fit to a logarithmic mass dependence using HV 
data above 3 x 10^^ h~^MQ and Virgo data below yields 

(c^) = (0.631 ± 0.001) (1- (0.023 ±0.002) In Af), (6) 
(c^) = (0.704 ± 0.001) (1 - (0.021 ±0.001) InM) (7) 

where M is the halo mass measure M200 expressed in 
units of lO^^h-^MQ. Although these expressions are 
good fits to the shape dependence over two orders of 
magnitude in mass, the limit c < 1 requires that the 
shape at masses approaching M200 < 10^ /i^^Mq must 
deviate this form. Simulations of smaller-scale structure 
will be needed to probe this regime. 

The rCDM model halos are more ellipsoidal and 
display somewhat weaker mass dependence than their 
ACDM counterparts. More vigorous growth of linear 
perturbations in the rCDM model drives a higher fre- 
quency of halo mergers in this model (Lacey & Cole 
1994), and this may explain why the rCDM halos have a 
mean c that is '--^ 0.03 smaller than the ACDM value. The 
velocity shapes of the tCDM model present a more puz- 
zling result, in that the logarithmic slope 0.012 ± 0.001 
is significantly shallower that the rCDM mass behavior. 
We suspect that our correction for Poisson noise may 
be inadequate for the velocities in this case (the mass 
resolution in the rCDM Virgo run is poorer than that 
of the ACDM Virgo run, so corrections are larger for 



TABLE 6 

Mass Dependence of Halo Minor Axis Ratio'' 



Model Component £15(0) a 

ACDM Position 0.631 ± 0.001 0.023 ± 0.002 

Velocity 0.704 ±0.001 0.021 ± 0.001 

tCDM Position 0.599 ±0.001 0.017 ± 0.001 

Velocity 0.692 ± 0.001 0.012 ± 0.001 

Trom z = halo samples. 

this model). The fiat behavior of velocity shape for the 
few most massive HV clusters also drives down the slope. 
The slopes of position and velocity shapes for the ACDM 
model are consistent; massive halos in this model have a 
fixed ratio of minor axis shapes c^/c^ = 0.896 ± 0.003. 

The behavior at z = motivates the following form for 
the behavior of the mean minor axis ratio at arbitrary 
redshift 

(c)(M,z) = Si5(z) (1-alnM). (8) 

Parameters at z = are listed in Table 101 

Turning to the redshift behavior, we first use the sky 
survey data to verify that the mass slope a does not de- 
pend on redshift. Binning clusters in two broad redshift 
intervals, from z = to 0.5 and 0.5 to 1, we fit the mass 
dependence of the mean minor axis ratio and find con- 
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Fig. 6. — The redshift dependence of cluster shapes, expressed 
in terms of the characteristic shape at lO-"^^ h~^MQ, derived from 
the combined HV sky survey outputs of the ACDM (circles) and 
rCDM (triangles) models. Filled and open symbols are velocity 
and position values, respectively. Lines are fits to equation j5J, 
with parameters given in Table 171 

sistency with the present-epoch slope. For example, the 
ACDM position minor axis ratio has a — 0.024 ± 0.003 
and 0.025±0.003 for the lower and higher redshift ranges, 
respectively, both of which arc consistent with the z~0 
value a = 0.023 ± 0.001. 

We characterize the redshift behavior of shape by fit- 
ting the mass-intercept at 10^^ h~^MQ to a power law in 
expansion factor 

ci5{z) = C15.0 (l + z)-". (9) 

In the sky sm'vey samples, each halo of mass M and 
minor axis shape c (bias corrected for Poisson noise) at 
redshift z contributes c/(l — alnM) to the estimate of 
Ci5(z). We use values of a measured at z = 0. 

Figure El shows the results derived from the combined 
sky survey samples (PO, NO, MS, and VS) of halos with 
M200 > 3 X lO^''/i"^M0, binned in Az = 0.1 intervals. A 
total of 44,122 (ACDM) and 19,813 (tCDM) halos are 
above the mass limit. Filled symbols in Figure are 
position while open show velocity shapes. Lines give fits 
to equation and the fit parameters are listed in Table 

HI 

The trend in redshift confirms the expectation that 
high redshift halos are slightly more ellipsoidal than their 
counterparts today. The redshift dependence is typically 
weak, £^^0.05 — 0.09, with the strongest trend exhibited 



TABLE 7 

Redshift Dependence of Halo Minor Axis Ratio" 



Model Component N^i £15,0 ^ 

ACDM Position 44122 0.633 ± 0.001 0.086 ± 0.004 

Velocity 0.701 ±0.001 0.053 ± 0.003 

tCDM Position 19813 0.600 ± 0.001 0.058 ± 0.007 

Velocity 0.691 ± 0.001 0.042 ± 0.006 

''From combined sky survey samples. 

by the position minor axis of the ACDM model. 

The present-epoch mean shapes at 10^^ H^^Mq are 
measured independently from the z = and sky survey 
datasets. The correspondence between £15(0) and C15.0 
to nearly three significant digits implies that our statisti- 
cal uncertainties are well understood and that there are 
no systematic differences at this level between the light- 
cone output used to generate the sky surveys and the 
more common output produced at fixed proper time. 

3.5. Comparisons to Previous Work 

Our shape results are generally consistent with previ- 
ous work, but different techniques for identifying halos, 
different shape measurement conventions, and variations 
in assumed cosmology complicate attempts at detailed 
comparison. 

The most recent work that is well aligned with our 
study is the N-body and particle hydrodynamics simula- 
tions of ACDM structure by Suwa et al. (2003). Their 
cosmological model has slightly higher normalization, 
(78 = 1, but is otherwise identical to that assumed in 
the HV and Virgo ACDM runs. They use a spheri- 
cal halo definition of somewhat larger radius (defined 
by an interior density contrast of 200 with respect to 
the mean, not critical, mass density) and find an aver- 
age minor axis ratio (c) =0.62 for clusters more massive 
than 2 x lO^'*/i~^M0. Their mass resolution is an or- 
der of magnitude improvement over the HV simulations, 
but their sample contains only 66 objects. Floor et al. 
(2003) use an annulus method to define projected shapes 
of halo from several hydrodynamical and N-body simu- 
lations. They find an average projected axis ratio of 0.77 
at z = 0, a value consistent with the mean intermediate 
axis ratio we find in three dimensions. 

Other published works employ an ellipsoidal region 
when defining halo shapes. Warren et al. (1992) ap- 
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plied this technique to ground-breaking, massively par- 
allel simulations and found a distribution of minor axis 
ratios for galactic-scale halos that peaked at c~0.6, with 
dispersion ^0.1. Thomas et al. (1998) analyze four Virgo 
simulations, including the two used here, and find mean 
minor axis ratio 0.50 with dispersion ~0.15 for a sample 
of 300 halos with mass limit 2 x 10^'* H^^Mq in low-f7,„ 
cosmologies. This shape result is significantly more ellip- 
soidal than our mean value of 0.64. The difference lies in 
the shape definition; we use moments of material within 
a sphere while they use moments of material defined us- 
ing a percolation algorithm on a set of particles whose 
local densities lie above a critical threshold. The bound- 
ary constraint of our method will tend toward rounder 
measures while the latter method, because of the direc- 
tional nature of the percolation process and the prun- 
ing in local density, will allow more strongly ellipsoidal 
values. The degree of difference between the two meth- 
ods can be large in extreme cases. For the same ACDM 
Virgo simulations, the most extreme position axis ratio 
measured by Thomas et al. (1998) is < 0.2, while the 
minimum value we derive is a factor two larger 0.4. Note 
that the effect in the mean is much smaller, ~ 20%. 

A more complete analysis involves measuring differen- 
tial shapes as a function of some scale parameter. War- 
ren et al. (1992) used an iterative scheme that began 
within spheres of fixed physical radius and found a high 
degree of correlation in direction and shape between 10 
and 40 Mpc in a sample of galactic halos. With much 
higher resolution simulations, Jing & Suto (2002) employ 
a sophisticated approach that first measures a local den- 
sity using a spherical kernel (as in SPH methods, Evrard 
1988; Hernquist & Katz 1989), then fits ellipsoids to par- 
ticles within some narrow range of local density. 

From 512'^-particle simulations of a ACDM cosmology, 
Jing & Suto (2002) measure the distribution of minor 
axis shapes at a density 2500pc (corresponding to a ra- 
dial scale of roughly r2oo/3) in a sample of 2494 halos 
more massive than 6 x lO^'^ /i^^Mq. They find mean 
(c) — 0.54 and dispersion 0.11 at z = 0. From analy- 
sis of twelve, high-resolution individual halo simulations, 
they find axial ratios that are rounder at lower densi- 
ties, with c~ (p/yOc)"'' ^^^- Using this relation to roughly 
scale to the radius used in this work, p/pc = 80 (equiv- 
alent to a mean interior density contrast of 200 for a 
p (X r~^-^ profile), results in (c) — 0.65 for their mass- 
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limited sample. However, our measurement of the mass 
dependence of shape implies that a second correction be 
made in order to infer their expectations at our mass 
hmit of 3 X 10^** /i^^Mq. Using equation (jSJ, the result 
is an expected mean (c) = 0.56. That the HV sam- 
ple mean of 0.64 is significantly larger may again simply 
reflect the different geometries being used in the two ap- 
proaches, but this hypothesis remains to tested. Note 
that the scaled Jing & Suto (2002) result is 0.05 larger 
than the mean derived by Thomas et al. (1998) for halos 
more massive than 2 x 10"'^'* h~^MQ. 

Clearly, there is not a unique measure of absolute halo 
shape, and future work is needed to more firmly establish 
the connections between different approaches to shape 
measurement. 

The trends with mass and redshift of the mean minor 
axis ratio presented in Figures |31 and are in qualitative 
agreement with Jing & Suto (2002). They find that halos 
of higher mass (at fixed epoch) and higher redshift (at 
fixed mass) tend to be more elongated. Their fit to the 
mass dependence is equivalent to a value a = 0.02, in 
good agreement with our finding of 0.023 ± 0.002. 

Regarding observed trends of cluster shape with red- 
shift, both Phonis (2002) using 903 APM clusters and 
Melott et al. (2001) using several optical and X-ray clus- 
ter samples find trends toward higher ellipticity at in- 
creasing redshift. However, Plionis et al. (2004) find a 
trend of shape with cluster size that is opposite that seen 
in simulations. From 1168 groups in the UZC-SSRS2 
galaxy group catalog, Plionis et al. find that poorer 
groups are more elongated than richer groups, with 85% 
of poor groups having b < 0.4. The discrepancy with the 
models may be due to biases in the optical group cata- 
logs or it may have a physical origin, such as galaxies not 
fairly sampling the dark matter in low mass systems. 

4. SHAPE AND LARGE-SCALE STRUCTURE 

For the linear initial density field, a connection between 
cluster shapes and large-scale structure was established 
by Bond, Kofman, & Pogosyan (1996), who showed that 
peaks separated by distances of order the mean inter- 
peak separation or smaller are likely to have strong con- 
necting filaments. By directing mergers occurring on op- 
posing ends, filaments serve as a source of alignment for 
halo shapes. Simulations show that this alignment per- 
sists into the non-linear regime (van Haarlem & van de 
Weygaert 1993). 
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In this section, we present mark correlations of shape, 
using two measures employed previously by Faltenbacher 
et al. (2003). We also examine whether clusters that are 
members of superclusters have shapes that differ from the 
general population. For the sake of brevity, we present 
results for the ACDM case only. 

4.1. Spatial Correlations of Shapes 

Several observational studies have probed the signif- 
icance and length scale of spatial shape correlations of 
galaxy clusters. Plionis (1994), using galaxy positions in 
637 Abell clusters, presents evidence for nearest neigh- 
bor alignment extending to separations 15 Mpc at ^ 
2.5a significance, with weaker alignment to 60 h^^ Mpc. 
Fuller, West, & Bridges (1999) used the brightest cluster 
galaxies (BCG) in poor MKW and AWM clusters and 
found significant alignments for BCG-cluster pairs and 
BCG-nearest cluster pairs. 

Simulations demonstrate that alignments are to be 
expected. Using a sample of several hundred clus- 
ters more massive than 1.8 x lO^^/i~^M0 derived from 
a 512^-particle ACDM simulation, Onuora & Thomas 

(2000) detect an alignment signal for pairs extending to 
30 h^^ Mpc, with the strongest signal for nearest neigh- 
bors. Like the position-velocity alignments of Figure O 
they note that the signal is persistent and changes little 
when only strongly elongated halos are used. 

Faltenbacher et al. (2003; hereafter F03) examine 3000 
clusters more massive than 1.4 x lO^'^/i^^M© from a 
ACDM simulation of a 500 h^^ Mpc region. They intro- 
duce the use of mark correlations to measure the behavior 
of alignment as a function of radial scale, and find pos- 
itive signal extending to 15 h^^ Mpc and 100 h^^ Mpc 
for two different alignment measures. Hatton and Ninin 

(2001) measure spatial alignments of halo pair angular 
momentums that extend to similar scales. 

Following F03, we quantify orientation alignment with 
two different measures. Consider an ensemble of clusters 
pairs, each with members i and j that have comoving 
spatial separation rij between r and r + dr. A basic 
alignment measure (termed A{r) by F03) uses the scalar 
product of major-axis directions d 

u{r) - {\d,-d,\) (10) 

where ( ) denotes an ensemble average over pairs of sep- 
aration r. 

A second, "filamentary" measure {J-{r) of F03) com- 
pares halo major-axis orientation to pair separation di- 




r (fi"' Mpc) 

Fig. 7. — Mark correlation of cluster orientations as measured 
by the mean excess basic alignment &u (open circles) and the mean 
excess filamentary alignment &w (filled circles) for the HV ACDM 
sample with A/200 > 3 X 10^'' /i^^Mq at 2 = 0. Solid lines show the 
fits to equations I12i and lliji . Filled squares show the sample's 
spatial correlation function in the regime where ^(r) > while 
open squares show \S,(r)\ in the anti-correlated regime. The dotted 
line is a rough fit of ^(r) = (r/22h~^ Mpc)~^-* within the range 
r ~ 20 - 6O/1-1 Mpc. 

rection 

w{r) - {\{\a^■h3\ + \dJ■h^))■ (11) 

For random halo orientations, both measures have ex- 
pectation value u{r) = w{r) = 0.5. We therefore use the 
deviations from this expectation 6u{r) =u{r) — 0.5 and 
Sw(r)=:w{r) — 0.5 as a measure of pairwise alignment. 

Figure [3 shows the ACDM position-space alignment 
signals at z = for the sample of 83,000 halos with masses 
-^-^200 > 3 X lO^*/i^^M0. Open circles show Su{r) and 
filled circles give Sw{r). Error bars are the uncertainty in 
the binned mean values. For comparison, squares show 
the two-point spatial correlation function ^(7') of the 
mass-limited sample with correlation length 22 Mpc 
(Colberg et al. 2000). Filled squares show positive spatial 
correlations while open squares show |^(r)| in the region 
of negative correlations (^(r) < 0). The inset shows u{r) 
and w{r) on a linear scale. 

The basic measure shows excess alignment out 
to scales ~ 30 Mpc. The filamentary statistic 
shows non-random halo orientations extending to scales 
200 h^^ Mpc, nearly ten times the correlation length and 
well into the anti-correlated regime of the two-point func- 
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tion. In both cases, the excess alignment signal is well 
fit by a power law at large radii that saturates at small 
r. The basic alignments follow 



5u{r) = 0.065 1 + 



-2.3 



while the filamentary alignment is well fit by 



5w{r) = 0.175 1 + 



-2.3 



(12) 



(13) 



2Ah-^ Mpc^ 

Our findings are generally consistent with those of F03, 
but our improved statistics and larger linear scale allow 
the first detailed fits to the effect. Although an analytic 
foundation for the specific form of equations 1)12(1 and 
1(13(1 is not yet in hand, we speculate that a solution may 
be found by applying the theoretical machinery describ- 
ing peaks in Gaussian random noise fields (Bardeen et 
al. 1986; Bond & Myers 1996). 

F03 note that the filamentary signal remains strong 
in projection, but their analysis is optimistic in that it 
assumes perfect knowledge of three-dimensional halo sep- 
arations as well as the projected three-dimensional halo 
shapes. A more appropriate treatment will require anal- 
ysis of galaxy cluster samples derived from mock catalogs 
(Kochanek & White 2003; Wechsler et al. 2004). 

4.2. Supercluster Members 

One might reasonably expect local cluster environment 
to play a role in determining halo shapes. In particular, 
since superclusters — groups of cluster-mass halos — 
tend to have strong filamentary morphology, one might 
anticipate that the formation dynamics of supercluster 
members may lead to a distribution of shapes that is 
biased toward more elongated systems. 

To address this question, we identify superclusters in 
the z = HV halo population by applying a percola- 
tion algorithm with linking length 23 Mpc, one-third 
the 69 h^^ Mpc mean intercluster spacing of the sample 
mass-limited at 3 x 10^"' H^^Mq. We further require that 
each supercluster have five or more halo members above 
the mass cutoff. This selects the 8% most strongly clus- 
tered halos. 

Figure IHl shows the distribution of minor axis ratios in 
position and velocity for the supercluster members and 
the general population. The axial ratios of the superclus- 
ter population have a nearly Gaussian distribution with 
means and dispersions presented in Table |S1 Comparing 
supercluster members to the general population (Table 
0)), we find no difference in mean values at the level of 0.01 
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Fig. 8. — Distributions of minor axis shapes for the 6683 halos 
of the ACDM supercluster population (histogram) compared with 
Gaussian fits to the general population (solid lines) at 2 = 0. A 
minimum mass of 3 X 10^"* h~^M.Q is employed. 



TABLE 8 

Supercluster Population Shapes'* 



Model 




c 


eg 


i, cr^ 


Position 


6683 


0.638 


0.092 


0.776 0.096 


Velocity 




0.704 


0.080 


0.818 0.088 



''Prom ACDM z = halo sample. 

in axial ratio; both have ((c), (b)) — (0.64,0.78) in posi- 
tion and (0.70, 0.82) in velocity. This finding lends sup- 
port to the picture in which the halo formation history 
is largely independent of large-scale environment (Bower 
1991; Sheth and Tormen 1999, 2004). 

5. CONCLUSIONS 

We use dark matter halos samples from billion-particle 
N-body simulations to measure statistical properties of 
halo shapes. The main ACDM samples consist of 83, 000 
halos a,t z = and 44, 000 halos from sky survey output 
with masses M200 > 3 x 10^'^ H^-^Mq. For each halo, a 
moment analysis of density and velocity structure within 
r2oo is used to define the principal axes of an equiva- 
lent ellipsoid. Higher resolution simulations and random 
realizations of isothermal spheres demonstrate that sys- 
tematic errors due to discreteness are less than 0.02 in 
mean axial ratio for the main HV samples. 

Massive halos have a Gaussian distribution of axis ra- 
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tios, with intrinsic dispersion ^ 0.08 and a mean minor 
axis ratio (c)(M, z)—ci5,o (1 — aln Af) (1 + 2:)^ that tends 
weakly toward rounder systems at lower mass and red- 
shift {a ^ 0.02 and e ~ 0.09 for position, see Tables 
and 13). 

Halos are rounder in velocity than in position space, 
a finding that likely reflects more efficient scattering in 
velocity space and the rounder nature of the gravitational 
potential compared to the density field. The principal 
axes in position and velocity are strongly correlated; 50% 
of halos have alignment angles smaller than 22 degrees. 
We also provide a Gaussian fit to the joint probability 
density of minor axis ratios in position and velocity. 

We investigate mark correlations of cluster shape using 
two statistics introduced by Faltenbacher et al. (2003). 
We measure significant excess alignment of halos ex- 
tending to 30 h^^ Mpc for the basic measure and to 
200 Mpc for the filamentary statistic. The latter 
is well fit as a function of scale by 6w{r) = 0.175[1 + 
(r/24/i~^ Mpc)]^^-^. The filamentary alignment should 
be detectable in large galaxy cluster surveys such as 
SDSS, but precise comparison between theory and ob- 



servation will require a careful study of the connections 
between clusters observed in redshift/color space and the 
underlying halo population. 

We also find that the distribution of halo shapes in su- 
percluster regions is indistinguishable from that of halos 
in general, showing that shape is largely independent of 
large-scale environment. 

The statistics presented here can be used to extend 
approximate methods for creating non-linear representa- 
tions of the matter distribution based on the halo model 
description (Scocciamarro and Sheth 2002) by incorpo- 
rating an ensemble of ellipsoidal halos with correlated 
orientations and aligned position and velocity ellipsoids. 
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